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Abstract 

We show how the introduction of a finite baryon density may trigger spontaneous parity violation 
in the hadronic phase of QCD. Since this involves strong interaction physics in an intermediate energy 
range we approximate QCD by a a model that retains the two lowest scalar and pseudoscalar multiplets. 
We propose a novel mechanism based on interplay between lightest and heavy meson states which cannot 
be realized solely in the Goldstone boson (pion) sector and thereby is unrelated to the one advocated 
by Migdal some time ago. Our approach is relevant for dense matter in an intermediate regime of few 
nuclear densities where quark percolation does not yet play a significant role. 

be" 

§ ' 1 Introduction 

Some time ago it was proved quite rigorously in pQ that parity, P, and vector flavor symmetry could not 
undergo spontaneous symmetry breaking in a vector like theory such as QCD. This is thus a well established 
result in strong interactions at zero chemical potential. Finite baryon density however results in a manifest 
breaking of CP invariance. The presence of a finite chemical potential leads to the presence of a constant 
imaginary zeroth-component of a vector field and the partition function of QCD is not anymore invariant 
under a CP transformation. The conditions under which the results of p] were proven (positivity of the 
measure) then do not hold anymore. 

The appearance of P violation for sufficiently large values of the chemical potential, i.e. at finite baryon 
density it is thus a logical possibility (conjectured by [2] in nuclear physics long ago). Can this possibility 
be realized in nature? Which would be its observable consequences should it occur? 

In order to answer the first question we could always appeal to lattice QCD for help. In fact, this possi- 
bility has been studied intensively for quite some time. However, finite density simulations are notoriously 
difficult (the fermion determinant is in general complex for a non-zero chemical potential fi) and one has 
to resort to more involved techniques such as determining the phase of the determinant separately, Taylor 
' expansions in [i or analytic continuation to imaginary chemical potential. On top of that Pauli blocking is at 

work for large values of //, etc. For a comprehensive review of the ongoing issues, see for instance [3]. Thus 
the lattice results for sufficiently large values of the baryon chemical potential (where the effect is expected 
to appear) are not known rigorously yelQ. 
Q\ ■ It is simpler to consider the (less relevant for physics) case of isospin chemical potential at zero baryon 

density [5] . In this case the fermion determinant is positive and spontaneous parity breaking is excluded by 
the Vafa-Witten theorem, although a pion condensate is in principle still possible. 

In this work we shall attempt to explore the interesting issue of P-parity breaking employing effective 
lagrangian techniques in the range of nuclear densities for which hadron phase persists and quark percolation 
does not occur yeo Our effective lagrangian is a generalized linear a model where we only include the lowest 
lying resonances, those that are expected to play a role in this issue. The use of effective Lagrangians is 
also crucial to answer the second question of interest, namely how would parity breaking originating from a 
finite baryon density eventually reflect in hadronic physics. 

Let us mention here several possible signatures of P-parity breaking, 
a) Decays of higher-mass meson resonances (radial excitations) into pions. Resonances do not have a definite 
parity and therefore the same resonance can decay both in two and three pions (in general into even and 
odd number of pions). 



< 
00 



Oh 



> 



o 



* On leave of absence from V.A. Fock Department of Theoretical Physics, St. Petersburg State University, Russia; 
andrianov@ecm.ub.es 
tespriu@ecm.ub.es 

1 The existence of a lattice parity breaking phase at strong coupling -the so-called Aoki phase- is well established and 
it is due to the fact that even at vanishing chemical potential the discretization of the fermion action leads to non-positive 
determinants [4|. This is clearly a lattice artifact and it is not the effect we are after. 

2 A qualitative estimate of the densities which enforce the overlapping of pions clouds, but not the cores, of nucleons indicates 
the necessity of including heavy meson states. 
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b) At the very point of the phase transition leading to parity breaking one has six massless pion-like states. 
After crossing the phase transition, in the parity broken phase, the massless charged pseudoscalar states 
remain as Goldstone bosons enhancing charged pion production, whereas the additional neutral pseudoscalar 
state becomes massive. 

c) Reinforcement of long-range correlations in the pseudoscalar channel and, correspondingly, additional 
isospin breaking effects in the pion decay constant and substantial modification of F v i for massless charged 
pions, giving an enhancement of electroweak decays. 

In the next section we shall introduce the model and see that there are some subtleties associated to the 
choice of the low-energy effective hadronic theory. Too simple models are not rich enough to explore all the 
different phases that the presence of manifest CP violation opens for us. We shall impose on the model the 
conditions for it not to lead to spontaneous breaking of parity for vanishing baryon chemical potential. In 
section [3] we shall introduce the finite chemical potential and see how it modifies the effective theory and 
the vacuum state. Next, in section [4] we shall understand how the masses and couplings of the particles 
are modified in such conditions. A specific model, that describes low-energy QCD, is studied in section [5] 
and it is seen that the values of the low-energy constants in QCD are compatible with the emergence of 
the phenomenon of spontaneous symmetry breaking at intermediate densities ranging from 3 to 5 times the 
normal nuclear density, approximately. The range of intermediate nuclear densities is of high interest as they 
may be reached in both compact stars [6l [7] and heavy-ion collisions [8] . 

There are some previous studies dealing with the problem of strong interactions at zero temperature and 
finite chemical potential: depending on a value of nuclear density, a variety of methods arc involved using 
meson- nucleon [21 H] or quark-meson [TU] Lagrangians and models of Nambu- Jona-Lasinio type [TTJ [T^] . 
One has also to mention an adjacent phenomenon of (C)P-parity breaking in meta-stable nuclear bubbles 
created in hot nuclear matter [13] . 

2 Generalized sigma model 

The simplest hadronic effective theory is the linear a- model of Gell-Mann and Levy [14] . which contains a 
multiplet of the lightest scalar a and pseudoscalar 7r a fields. Spontaneous chiral symmetry breaking emerges 
due to a non-zero value for (a) ~ (qq) / F% . Current algebra techniques indicate that in order to relate this 
model to QCD one has to choose a real condensate for the scalar density, with its sign opposite to current 
quark masses, and avoid any parity breaking due to a v.e.v. of the pseudoscalar density. The introduction 
of a chemical potential does not change the phase of the condensate and therefore does not generate any 
parity breaking. This is just fine because in normal conditions P breaking is impossible in QCD. 

However, if two different scalar fields condense with a relative phase between the two v.e.v. 's the oppor- 
tunity of spontaneous parity breaking may arise. 

Let us consider a model with two multiplets of scalar/pseudoscalar fields 

//, ^1 • J = 1,2; HjH} = (a] + (ir?f)I, (1) 

with itj = TTjT a with r a being a set of Pauli matrices. We shall deal with a scalar system globally symmetric 
in respect to SU(2)l x SU(2)r rotations and work in the exact chiral limit, with vanishing current quark 
masses. We should think of these two chiral multiplets as representing the two lowest-lying radial state^j for 
a given J PC . Of course one could add more multiplets, representing higher radial and spin excitations, to 
obtain a better description of QCD, but the present model already possesses all the necessary ingredients to 
study spontaneous parity breaking (SPB). Inclusion of higher- mass states is required at substantially larger 
densities when typical distances between baryons are shrinking considerably and meson excitations with 
Compton wave lengths much shorter than a pion one play an important role. 

Let us define the effective potential of this generalized a model. First we write the most general Hcrmitian 

3 Some previous attempts to give the Lagrangian description for two multiplets of scalar and pseudoscalar mesons have been 
undertaken in |15j (the oldest one ) and in |16| (the most recent one). However we have been rather inspired by our previous 
works on extended quark models |17II18I . 
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potential at zero \x compartible with SU(2)l x SU(2) 



r 



V eS = itr J - ]T HtA jk H k + Xi{h\Hi) 2 + X 2 (h\h 2 ) 2 + \ 3 h\HiH\h 2 

+ ^X 4 (hIh 2 h\h 2 + HlHxHiHx) + \x 5 {h\h 2 + h\Hi)H\Hi 

+ \x 6 (h\h 2 + HlH^HlH^ + 0(^f), (2) 

with 9 real constants Ajk, Xa ■ QCD bosonization rules indicate that Ajk ~ A^ ~ Ar c . The neglected terms 
will be suppressed by inverse power of the chiral symmetry breaking (CSB) scale A ~ 1.2GeV. If we assume 
the v.e.v. of Hj to be of the order of the constituent mass 0.2 0.3 GeV, it is reasonable to neglect these 
terms. In section [5] we shall estimate the values of the constants A,4 approximately corresponding to QCD. 

After specifying the v.e.v. (Hi) = (ci), one can use the global invariance of the model to factor out the 
Goldstone boson fields with the chiral parameterization 

Hi(x) = ai(x)U(x) = ai(x)Z 2 (x); H 2 (x) = Z(x)(a 2 (x) + Ht 2 {x))S(x). (3) 

This kind of parameterization preserves the parities of a 2 (x) and tt 2 to be even and odd respectively (in the 
absence of SPB) and realizes manifestly the masslessness of Goldstone bosons. 

Let us now investigate the hypothetical appearance of a non-zero v.e.v. of pseudoscalar fields. In order 
not to break the charge conservation, we must expect, if at all, only a neutral condensate represented by a 
solution with tt 2 — 5 a0 p. The conditions to have an extremum are derived from the first variation of the 
effective potential ^ after substitution of Eq. ^ , 

2(AuCTi + Ai 2 ct 2 ) = 4Aicri i + 3A 5 cti er 2 + 2(A 3 + Xi)aio 2 2 + X & a\ + p 2 (2(X 3 - X 4 )ai + A 6 cr 2 ) , 

2(Ai 2 cri + A 22 cr 2 ) = Ascr? + 2(A 3 + A 4 )ct?(t 2 + 3A 6 ctict^ +4A 2 cr^ +p 2 ^A 6 cri +4A 2< t 2 ), 

= 2n a 2 ( - A 22 + (As - Xi)a 2 + X 6 a x a 2 + 2X 2 a 2 + 2A 2( o 2 ) . (4) 

To avoid spontaneous parity breaking in normal vacuum phase of QCD, it is necessary and sufficient to 
impose 

(A 3 - A 4 )(7 2 + A 6 (7i(7 2 + 2A 2( t 2 > A 22 . (5) 

Since QCD in normal conditions does not lead to parity breaking, the low-energy model must necessarily 
fulfill©. 

A necessary condition to have a minimum for non-zero <Tj (for vanishing p) can be derived from the 
condition to get a local maximum (or at least a saddle point) for zero <jj. The sufficient conditions follow 
from the positivity of the second variation for a non-trivial solution of the two first equations J3J at p = 0. 



3 Finite chemical potential 



We shall assume that the scalars under consideration are generated in the quark sector of QCD. The baryon 
chemical potential is transmitted to the meson sector via a quark-meson coupling. Without loss of generality 
we can choose the collective field having local coupling to quarks as Hi; this actually defines the chiral 
multiplet Hi. The set of coupling constants in ^ is sufficient to support this choice as well as to fix the 
Yukawa coupling constant to unity. Thus finite density is transmitted to the boson sector via 



AC = -(q R Hiq L + q L H\q R ), 
where qL,R are assumed to be constituent quarks. Then the one-loop contribution to V e s is 



(0) 



avm 



\Hi\) 



p\Hi\ 2 ^/p 2 ~\Hi\ 



|V-|ff l | a ) 3/a 



#i In 



, p+^p 2 -\Hi P 



\Hi\ 



(i + o( 



d \_Hii 

A 2 ' A 2 



Af : 



N r N 



4tt 2 



(7) 



3 



where p is the chemical potential. The higher-order contributions of chiral expansion in 1/A 2 are not 
considered. This effective potential is normalized to reproduce the baryon density for quark matter 



PB = -\d,AVM = ^f P % = ^( M 2 - \{H^f/\ (8) 



where the quark Fermi momentum is pf — \/ P 2 ~ \(Hi)\ 2 . Normal nuclear density is ps — 0.17 fm~ 3 ~ (1.8 
fm) that corresponds to the average distance 1.8 fm between nucleons in nuclear matter. 

After the introduction of p the conditions for a minimum of the effective potential are modified. In 
particular, the first equation in Q takes the form 



2(A n oi + A12CT2) = 4Ai<t 3 + 3A 5 cr 2 cr 2 + 2(A 3 + A 4 )ctict^ + A 6 cr| + p 2 (2(A 3 - A 4 )cti + A 6 cr 2 
+2MQ(p - <7i ) 



pui \ p 2 — o\ — o\ In 



p+ yfp 2 - a\ 



01 



(9) 



and other ones are unchanged. 

Let us first focus on the regime of small chemical potentials. For small values of p, we know that the 
value of the odd parity condensate p is zero. The possibility of SPB is controlled by the inequality in 
order to approach a SPB phase transition we have to diminish the l.h.s. of inequality ([5]) and therefore we 
need to have (assuming that the inequality indeed holds at p = 0) 



(A3 - A 4 )(T 2 + A 6 (Tl(T2 + 2A 2 CT2 < 0. (10) 



Let us now leave the case p ~ (01) and examine the possible existence of a critical point where the strict 
inequality §5§ does not hold and instead for p > p cr n 

(A 3 - A 4 )(J 2 + A 6( ji(72 + 2A 2 (erf + p 2 ) = A 22 . (11) 

After substituting A22 from (jTTJ) into the second equation in (HJ one finds that 

A 5 CT 2 + 4A 4 cricr 2 + Aef^ + P 2 ) = 2A12, (12) 



where we have taken into account that <ti ^ 0. Together with (jTYJ) this completely fixes the v.e.v.'s of the 
scalar fields 0^2 . If A2 = and/or Xq = the corresponding Eq. (jTTJ) or (JT3J) firmly fix the relation between 
and (72- Otherwise if A2A6 ^ these two equations still allow to get rid of the v.e.v. of pseudoscalar field. 
Thus in the P-breaking phase the relation between the two scalar v.e.v's is completely determined. Using 
Eqs. (J9j) and (fl2|) one can easily eliminate the variables p and 172 and get for numerical calculations the 
equation for o\ solely. 

Let us now try to determine the critical value of the chemical potential, namely the value where p(p c ) = 0, 
but Eqs. fTTT) . ([T2l hold. Combining these two equations, 

(4A 2 Ai2 - A 6 A 22 )r 2 + (2A 6 Ai 2 - 4A 4 A 22 )r + 2(A 3 - A 4 )A 12 - A 5 A 22 =0; r = — . (13) 

01 

In order for a SPB phase to exist this equation has to possess real solutions. If 4A2A12 — A6A22 ^ the 
SPB phase is bounded by two critical points corresponding to second order transitions. If, on the contrary, 
4A 2 A 4 2 — A@A 2 2 = there is only one solution corresponding to a second order transition, but there may 
exist other solutions that fall beyond the accuracy of our low energy model (which becomes inappropriate 
for small values of 01). See section [S] 



4 The physical spectrum in the SPB phase 

Once a condensate for 7r° appears spontaneously the vector SU(2) symmetry is broken to U(l) and two 
charged tt' mesons are expected to possess zero masses. In this case the matrix of second variation 
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reads 



■y(2)<r 
! 11 

y(2)a 
v 12 



(2)c 



4Aio-f + 2A 5 f7iCT2 + 2A40-2 - 27Vc7i In 

2A5O-1 + 4A 3 f7lCT2 + 2A 6 (T2, 

4(A 3 - A 4 )cri +2A 6 cr 2 )p, 



+ JpT~o{ 



CTl 



:1 



(2)<r 



2A 4 erf + 2A 6 cricr 2 + 4A 2 ct 2 , 
V$™ = (2A 6 a 1 +8A 2( r 2 )p, 



^00 



(2)tt 



4A 2( o 2 



1 



(2)tt 



0. 



(14) 



where the r.h.s. are evaluated with the help of Eqs. lfTTj) and (fl2"|) . We notice that convexity around this 
minimum implies that all diagonal elements are non-negative Vjp a > 0. This gives positive masses for two 
scalar and one pseudoscalar mesons, whereas the triplet of pions and charged doublet of ir' mesons remain 
massless. Of course, the mass spectrum can be obtained quantitatively after kinetic terms are normalized. 

Once we have fixed the interaction to quark matter we are not free in the choice of the kinetic term for 
scalar fields. Namely one cannot rotate two fields and rescale the field Hi without changes in the chemical 
potential driver (0. However the rescaling of the field if 2 is possible at the expense of an appropriate 
redefinitions of other coupling constants and this freedom can be used to fix one of the constants which 
appear in the kinetic term. Thus we take the general kinetic term symmetric under SU(2) l X SU(2) r global 
rotations to be 



1 



(15) 



j,k=i 



After selecting out the v.e.v. (Hi) = (ci) = ui one can separate the bare Goldstone boson action with 
the chiral parameterization Q. Let us explore the kinetic part quadratic in fields. We expand U = 



1 + m/F + • • ■ , i = 1 
part looks as follows 



i7r/2i 7 o + ■ • ■ and use the v.e.v. 's ctj = Oj + Sj 7r = r^p + II. Then the quadratic 



.(2) 



j,k=l 

tX a * 



1 , 
- 2 A 22 



(16) 



which shows the mixture between light and heavy pseudoscalar states and, in the SPB phase, also between 
scalar and pseudoscalar states. 
Let us define 



1 



3,fc=l 



3=1 



In the symmetric phase p = 0, (|16j) can be diagonalized, 



r( 2 ) — -F) n a r)!- l TT a 



1 



(A 22 -c 2 )^n a ^ir, a 22 -c 2 = 



cffdctA 



F 2 



>o, 



(17) 



(18) 



wherefrom, taking into account the matrix of second variations of the effective potential l[2|). one finds the 
masses of the heavy pion triplet 



m 1 



-A 22 + (A 3 - A 4 )(ai) 2 + A 6 aia 2 + 2A 2 (a 2 ) 2 



A 



22 



c 2 



(19) 



In the SPB phase the situation is more involved: pseudoscalar states mix with scalar ones. In particular, 
diagonalization is different for neutral and charged pions because the vector isospin symmetry is broken: 
SU(2)v — * U(l). The SPB induces mixing of both massless and heavy neutral pions with scalars. The 
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(partially) diagonalized kinetic term has the following form 

4i = d^d^ + + if^!w«W + (A 22 - c^n^iF 

+ 2 X, ?fTw ^ fc " fTw C9p 4- 329 " (20) 

— i j — i 

We see that even in the massless pion sector the isospin breaking SU(2)v — » U(l) occurs: neutral pions 
become less stable with a larger decay constant. Another observation is that in the charged meson sector 
the relationship between massless 7r and II remain the same as in the symmetric phase. 

5 P-violation in models with discrete Z 2 x Z 2 symmetry. 

As a relevant example we now examine models with residual discrete chiral symmetry (after the breaking 
SU(2) L x SU{2) R -> SU(2) V ) under independent reflections Hi -> -Hi and/or H 2 -> -H 2 ■ Then 
A5 = Ag = 0, A12 = 0, A12 = 0, but A4 ^ 0. One can always fix A\ = A 2 redefining the other parameters. 

Let us now see how the general relations in the previous section are realized in this model. The analysis 
of Eqs. (2]), © and ([TO]) as well as the positivity of the second variation matrix leads to conclusion that 
in the symmetric phase the only solution compatible with the very possibility of P-breaking is a 2 = . In 
the SPB phase for these models the constraint 1)1 2p 2A 4 eri<72 = A12 = has also a unique solution a 2 = 0. 
Therefore <j 2 = everywhere. As for ai we get for fi = 

2 A n 
ffl= 2V 

The condition for extremum when \x 7^ reads 



An - 2\\o\ - (A 3 - A 4 )p 2 = Af@(n - a t ) 



H\ hi 2 -o{-ai In 



2, vV - O- 



(21) 



0"1 

and the condition for this extremum to correspond to a SPB phase is now 

(A3-A 4 )a?+2A2P 2 = A 2 2. (22) 
The second variation matrix for 172 = in both phases reads 

)vg )a = 2A n - 2(A 3 - A 4 )p 2 - Wft J - u\ > 0, ^ 4A 1( x 2 , 



2 

V[T = 0, iy 2 ( 2 2)CT = -A 22 + (A 3 + A 4 )a 2 + 2A 2 p 2 S ™ 2A 4 a 2 > 0, 
= 4(A 3 -A 4 )a lP , F 2 ( 2 ^ = 0, 



iV^ 2)w = -A 22 + (A 3 - A 4 )a 2 + 6A 2( o 2 ^5 4A 2( o 2 > 0, 



-Vg^ = -A22 + (A 3 - A 4 )a 2 + 2A 2P 2 S -^5 0. (23) 



As this matrix is positive definite one derives the following requirements for the low-energy hadronic model 
to make sense 

Ai>0, A 2 >0, A 4 >0, An >fj, (A 3 ± A 4 )A U > 2AiA 22) (24) 

where the last inequality is obtained from the positivity of Vjjf , ^00 a t zero /i . 

As it follows from (113p there is now (at most) one critical point (within our approximations) where 



G 



For this solution to correspond to a real value for o\ and therefore be acceptable, the r.h.s. of (125j) has to 
be positive. In addition we need Eq. (f2"T|) (after substituting p = 0) to have a solution for /j, in the range 
f 1 < M < °°7 



(A 3 - A 4 )An - 2AiA 22 
(A3 - A 4 ) 



2 2 2 1 M + vV 
A*V M - o-{ - ct x ln- 



> 0. (26) 



0"1 

Then after comparing inequalities in (f2"4"]) , ([23]) and (|26[) one concludes that 

A 3 > A 4 > 0, A 22 > 0. (27) 

Thus the SPB phase should exist for all positive constants \ } > 0, Ajj > if A3 > A4. The critical value of 
chemical potential can to be calculated from Eq. (|26[) . 

Now let us estimate the typical scales of P-breaking from meson spectroscopy. Since we have defined 
A\ = Ai we can find from Eq. p3|) the ratios of masses. Just to get a feeling of the possible scales involved, 
let us make a tentative choice 

A n = - = A 22 , F Q = lOOMeV, 1 = 300MeV = 3F , 
9 

according to the relation Fq = A\o\ , and use the units Fq further on. As well 

772 7t = 0, 77i Sl = 0.7GeV = 7F a , m n = 1.3GcV = 13F , to S2 = 1.5GeV = 15F (28) 

in a fair agreement with particle phenomenology [19]. Then from the definitions of masses 

2 2A n 4Ai^ 2 2 2A 4 a? 
™ Sl = = —1 — ; ™s 2 - mn = —. — , (29) 

one finds 

An ^ 2.7F 2 , Ai ~ 0.15, A 4 ~ 0.35. 

Taking cri tCr it — l-Si^o and the previously estimated value Ai = 0.15 one finds that SPB occurs at 
Pf — 3.9Fo = lAApF,nudear which corresponds to dense nuclear matter with QB.crit ~ 0.5fm^ 3 ~ 3g B . 

nuclear • 

The phase transition occurs at /i c ~ 4.3Fo > o\. From the definition of (|25|) and of the mass of II one finds 
A 3 ~ 3.6 and A 22 ~ 11F 2 . Thus we see that the possibility of SPB emerges naturally for reasonable values 
of the meson physics parameters and low-energy constants. At this critical point the masses of scalar mesons 
are mjj ~ 1.7Fo, ms 2 ~ 4.5Fp . 



6 Conclusions 

Let us summarize here our main findings. Parity breaking seems to be quite a realistic possibility in nuclear 
matter at moderate densities. We have arrived at this conclusion by using an effective lagrangian for low- 
energy QCD that retains the two lowest lying states in the scalar and pseudoscalar sectors. We include a 
chemical potential for the quarks that corresponds to a finite density of baryons and investigate the pattern 
of symmetry breaking in its presence. We have found the necessary and sufficient conditions for a phase 
where parity is spontaneously broken to exist. In general this phase is bound and it extends across a range of 
chemical potentials that correspond to nuclear densities where more exotic phenomena such as color-flavor 
locking or color superconductivity may occur. 

Salient characteristics of this phase would be the spontaneous breaking of the vector isospin symmetry 
SU(2)v down to £7(1) and the generation two additional massless charged pseudoscalar mesons. We also 
find a strong mixing between scalar and pseudoscalar states that translate spontaneous parity breaking into 
meson decays. The mass eigenstates will decay both in odd and even number of pions simultaneously Isospin 
breaking can also be visible in decay constants. 

We think that our conclusions are drawn in a region of parameters where our effective lagrangian is 
applicable and, while obviously we cannot claim high accuracy in our predictions, we are confident that the 
existence of this novel phase is not an spurious consequence of our approach but a rather robust prediction. 
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It would surely be interesting to investigate how this new phenomenon could possibly influence the equation 
of state of neutron stars (the density of such objects seems to be about right for it). 

Lattice methods could shed some light on this issue and confirm or falsify the existence of this interesting 
phase in dense nuclear matter. One could probably use the expansions for small values of fj, at finite 
temperatures to check some of our expressions. For this matching the natural approach is a hot hadron gas 
|20j . Conversely, it would be possible to extend our techniques to the case of isospin chemical potential [5] . 

We are grateful to our collaborators S. Afonin and V. Andrianov for checking the calculations and 
making very useful remarks . This work was supported by research grants FPA2007-66665, 2005SGR00564, 
2007PIV10046. It is also supported by the Consolider-Ingenio 2010 Program CPAN (CSD2007- 00042). We 
acknowledge the partial support of the EU RTN networks FLAVIANET and ENRAGE and the Program 
RNP2. 1.1. 1112. 
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